<asm«m> k<s. 

•%!mS1£#l!im$.-fa77MZ08SmT-t>T7) ~ 75 

•#l?eifflfied*aS'7 , o^7^(BASIC) --- 77 

•^eftsawaoisyD^ACBAsio 79 

• 2 * • 3 *Zm£&&V®<t7 u? ? j, (BASIC) 95 

ma&fcwtwm-rar? a <basic> i 01 

•S»iaK»¥«Ul«rDjr?A(8085fflT*>T5) 107 

<e*»iB««Hi/aia»> $6$ 

•b--?«asrB^7A(BASic) • 113 

•«S«raB8»»t«VD /Ji (BASIC) U9 

<f-t > a >m> #7* 

••f <)a.— > a ^Jta7"Dr7A(BASIC) 129 

<8«7-uia:««> sjs 

• Sande-TukeyffiCJ:{FFTrD^7A(BASIC) 154 

• Sande-Tukey mtztl FFT7*d (8086JH7-fe >T5) • 15 6 

• tr •* h ae/sin/cos ?• -T^ar d *r ? a (basic) iei 

• FFT»fifrD77^ffll(BASIC) 163 

<**i>>ne-as> »9« 

• MEMyD^7Agj(BASIC) 174 

•to«'<7>-^jfcgffly D j'7A(BASic) 189 

• v-syisytxiktzi 4tt»fl-i8yD^5A(BASIC) 19 3 

• DFPteci4fcj&JM|7 , D#'?A(BASIC> 198 

• Gauss-Newton ffiC J: 4i£0fl«7'o y* 7 A (BASIC) 205 

6 ft* 2 ft^ro 7 5 A (BASIC) - .... n6 

•ia*ffiw#frD77A(BAsic) 222 

ffi: BASIC II, rc-980i«rtr>'V:JWSii>( ( »T>*i. 
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»ff<cirT-WSvS*iT(..« -/tj->Ba- atBRTSil. »5H 

«. «a«2, ffi*m«*MKHo5a*, aieattr, stztttBifiizism-ce*. 

11.1 «£gjs*r — m&tma 

«*<EKfc-*K«-(tta. fcffira)T-|ii<, i/.t*>o-MfiS-(t£«, fflB5)ft-3ife 
If, iii^*"»«"tjat*'-egs. oil), *<»(g«raiji*.t, Sffltii><-3*> 
«>«a*t>-5i>*lfT*s< * i*>, *M,m1#at*«fft4;:ifc:.t,T, £ttroffi«S* 
fflSLi-j it4rar*4. 

#il4rW* A If T A i 7 . £ i £.*tl?tUSm<Z&*tlt:B(l80 1 8085. * LT Z 80<* 4 



11.2 ftBgggwKB iu 

4f*«>T)ofta«s***o*<i'*r»a«*t. *fe»>6*ti..?ii«>ftatt«»att. 

■ SB - * 7 -o<0*S5».# JEgJW(Multi7ariate Analysis) 1 1 ti h 4 b o>V f, 4 . 

tse<oaa^«oBWcJ: ->x. ®mis&m. ©±b&hh«s. ©was8&& ® 
11.2 fi0*l£*fCDfi58 

*a<¥ffit*'*4S', aiasMrosiscH+i. fi®»#*re.t9Sfeft./i:'Hffi-ii. 

8>SS • ft* • Sft^W'jrttBSKIBt 4 S 4 ?4 <f mt*L 4 t OT» 9. S felcli. 
-»mSK*feMl*«3f|ffflfT'M-r4raict«mTS4. CiTIJ, &£T-t U-c? 

»ra*ttft3fc<>rax-t* ^*«i«i.-<* mxtohhmt-<l> i- 



Ell. I ffi£-tt*-<? KkBt^rMl 
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atftTi/ »dc-c* hn- c =(c,.<i. ",c)'-ca-r. fccu < wests-*-. ;:tis# 

= Sr+n 



<11-U 



Q-\*-8lV (11-2) 



-52.= o. -52- = 0 -, -2fi- = o 
at, ai , • az„ u 



(11-3) 



tWittfiv-. < to 4> h.T i- 4 i i fc 

e = (S'S)-'S'* (11-4) 

mm. m&m&m smth. 



U.3 KBB»t»ayoy9Z,tB«gq 115 
8 = {S'GS)-'S'G* (US 

*c tft*7 tEfcre 4 {. <or*h.«, c <ott*s?« *> a i/jr.(x,iia2E6-&A-< 

? t-zl-r <Oi*Bra5jS}>fcL. (iOjIflit-itOttWIlK 

S 0 ( i = 1 . Jf >/>ft**ftO 6 t T(ll-2)*<0«+ 2 ftRKe i !)» < . t* L, 

SJmiiCtcloT (11-4) Soii c«5$o»ra£»6# (ltli?»f, 
tUtli*SK«SK&4. i(lt»<7*f)Xi.KttK<-3*ti4i, t -> fc fc i < J8 
K(.tl4Wli. #a»ttatt(-*»taa)T*4. C»?&e-»<T<><fclA»||tq!UX 
K3)S»fflStl.fcv>. 

11.3 wmmti^yav^LtMmmi 



► BASIC Vay"}U 

affit4. ift.l±. S±aFM-8fflOFBASICTS*KTv>4*'. i <=> -i/t^Otf I 

*ffi)iri)>#(e*ai*m)tffimt4.Hii.2tc*flftJt£osa+2*ao7 , D^7A 

HU.2«37o77.tcm>T, »S-3SS4> 2 «ffit*nt4/w-fl' (1510- 2420}?) 
6ATffc«-5ftv». 

RKC3IJ)Ma*»JtlcJ:5«fi»«ff'-^lc«LTH11.2<oyDj'?AtHv.-C(S'&a 



Men Negative 

LMtt &qu«rw NatlKKl 



looo • 

1010 ' SSIt«t**llS***t(Stt»S«S *****••• ** 
102O -t 
1030 
two '* 
1030 *« 

1060 'ttlt*S*tSI«ttttt»«ttttt»SX»tfltl*S 

1070 * 
1080 ' 
1090 ' 
1 too ' 

1 1 lO 0EF1MT J-N 

1120 DEF FNZHt>-100-liE*o-T*(-2>- 
1130 ■ 

JI40 'i 1 

IISO ' I nai n Rout in* I 

1160 ' ' > 

1170 • 

iibo wiirm 00,23 

1190 C06UB 1330 

1200 Din Z(N3>,RfHC,NC>,BCNC>,X<NC> 
1210 OJM R14NC,NC>,B1<NC),XMNC) 
1220 Dili IZ<KCl,IP(NC),JZ(NC>,MtNC> 
1230 TIHE»--00i00i00- 
1240 80SUB 1SS0 
1250 TOTIHE* 
1260 GOSU9 2400 
1270 STOP 
12BO ' 

1290 * i ■ , 

1300 ' I Date. Input Rout in- I 
1310 * « - ,■ 
I320.* 
1330 CLB 

1340 PftlMT SPC(A) "**• NMLB tt*», PRINT 
1350 INPl/T'Nuaticr of Component "jNC 
13AO IW>UT*N(Mb«r o4 Sanple Point • j MS 
1370 Dill 80fNC,NSl, D<NS} t A*4NC> 
1380 FOR 1-1 TO NC 

1390 INPUT "File Name of Standerd " l A* ( 1 > 
1400 NEXT 

1410 INPUT-Flle Nan* of Gaapla -lAA* 

1420 FOR J=l TO NC 

1430 OPeN"I-,«l t -j,-+A*<J) 

1440 FOR 1-1 TO NSiINPl/T*l,6D4J,I>iNEXT 

14SO CLOSEel 

1440 NEXT 

1470 OPEN"l- t «l f -l l "*AAe 

14B0 FOR J»I TO M3j1NPIIT*I1,D(JI JNEIT 

1490 CLOSES 1 

13O0 RETURN 

310 • 
I320 



j fli*A-<7 * 7 <f Jt* A f] 



1530 ' I Non Negative 



Least Squares Method 

Execution Routine 



1340 : 
1330 ' 1 

1S40 ' < 

1370 * 
1SSO FOR 1 = 1 TO N6 
1370 Z(I)-FNZ(D<I>) 
14O0 MEXT 

1610 FOR 0-1 TO MC i FOR K"J TO NC 



1620 Rfi-O 

1630 FOR 1*1 TO NS 

1*40 RR-nR*CSU,I)*&D<K,I>/Z<I> 

1630 NEXT 

1660 R<J.K>-AR; R(K;j)=RR 

1670 NEXTiNEXT 

1680 FOR J»l TO NC 

1690 BB-O 

17O0 FOR 1-1 TO NS 

171© BB-BB+SOtJ,I)*0(l)/2(lt 

1720 NEXT 

1730 B<J>-B8 

1740 NEXT 

17S0 N2-NCiNPeO 

1760 FOR 1-1 TO MC 

1770 X<l>aOilZfI)»I. 

i78o we XT 

1790 FOR 1*1 TO NC 
1SO0 RX-0 

1810 FOR J=l TO NC 
1B20 RXHtXHHl, J)*X1J> 
1830 NEXT 
1940 W(I>-B(I)-RX 
1030 NEXT 

10-60 IF NZ-O THEN 2420 
1670 UMAX«0lMc0 
1BS0 FOR J-l TO HZ 
IB TO lloIZtl) 

I900 IF UIIDMMIII THEN WHAX-UM 1 ) »HJ=I 
1910 NEXT 

1920 IF HI-0 THEN 2420 

1930 NF*NP+liIPtKP)-lZ(MI> 

1940 lZ(ni)-JZtNZ) iNZ-NZ-1 

1930 FDR 1-1 TO NP 

1960 IloZPtJ) 

1970 FOR J-l TO NP 

I960 Jl-lP(J) 

1990 R](1,J)4< II, Ji> 
2000 NEXT 
2010 Bl ri)-S(It) 
2020 NEXT 

2030 IF NP-1 THEN XI (1)=BH1> /Rl <i f 11 iGOTO 2190 

2040 FOR 1-1 TO NP-liFOR 0-1*1 TO NP 

2030 GC-R1 40, I>/R| <i, ]) 

2060 FOR K-Z+i TO NP 

2070 RlfJ,K)-Rl(J,IO-6CtRl(X.K> 

2060 NEXT 

209O B1<J)-81(J)-SC*BMI) 
21 OO NEXTiNEXT 
2110 X1<NP1«B1<KP)/RWNP,NP> 
2120 FOR 1-NP-l TO 1 8TEP -1 
2ISO srt-O 

21 40 FOR 0-1*1 TO NP 
2130 EM=Sn+RKI,J)»Xl (J) 
2160 NEXT 

2170 X1<I)-<B1(I>-SM)/R1<I,I) 
21B0 NEXT 
21 TO J*Q 

2200 FOR 1-1 TO NP 

22IO IF XMIX-O THEN J-J+l 1 JZ (J » -I 
2220 NEXT 

2230 JF J<X> THEN 2260 
2240 FOR I »l TO NP 



IS But* 

2230 I<IP<1>>'IMI> 

2240 NEXT 

2270 B0TO 1790 

22BO ALKItllleO 

22VO FOR 1-1 TO J 

2200 H=J7(I)iI2=IP(H> 

2310 AL1-X<I2>/(X(12>-XKI1>> 

2520 IF ftLKAL THEN ftL-ALl i n I- 1 1 

2310 NEXT 

2M0 FDR 1=1 TO NP 
233© 11-IPU) 

2360 X<U)-(l-AL)4X<llI*Al.aXl«l> 

2370 NEXT 

23BO XfIP<HIJ>^> 

2390 NZ*N7«11IZ(NZ]clP(nl> 

24O0 IPCM>-JPlNPIiNP-NP-l 

2410 GOTO 1790 

2420 RETURN 

2430 ' 

2440 * i 1 

2430 ' I Output Routine I 

246© ' ' 1 

2470 ' 

24BO PRINTfpRINTaPftINT -Non Nwqatlva LMtt Squwas C E~2= 

•|E»j" 3-iPRIMT 
2490 PRINT EPCC71-0ATA i -|AA4|PR1NT 

2300 PRINT SPC<6>|iF0R 1-1 TO NCiPfllNT SPC<7» | If iNEXTt PRINT 
. TINE" 

2310 PRINT SPCI7)itF0R 1-1 TO NCt PRINT USINS-fttpt*. «•••■> X 

till (NEXT I PRINT ■ - 1 T« 
2320 OPEN"0-,«l,"LPTOi " 

2530 PRINT! 1. "Woo Nvgttlw LsaEt. Squares I E**2-"|E«|" 3-iPRINTal 
2340 PRINTI1,EPCI7I-0ATA ■ •|«A«iPRINT4I 

25SO PRINTal,SPC(6]|iFOR I«l TO HCi PRINTt 1 f SPC (T> 1 1 $ iMEXTtPRJNT»I , - 
TINE ■ 

25AO PRINTIl.SPC(7>|:FOR 1=1 TO NClPRINTIl ,U31KB-tl>»a. X II ) , 

iNEXTiPRINTU, - -|T« 
2570 RETURN 
2SSO * 
2S90 END 



tSStffoii. CWSS«r«, 3«*>eK(A : J"Hi'-it\'-sV. B :7n^9w/-n, • 

BU.4(«),fl>)lc-omig^feX„X,«i«(Sx-i# in-- f-*(300nm-700nm, .nm 
.Tfcl01j&)£zRT. X,li3a<O6g5t0.6:0.3:0.2Wit¥TiiiiiT*9, X.ljj'f-fl'-* 
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at* -< >9 - -3 * - X IEEE-488 J 
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(ftftUV-lBO) 



(X±Tft,FM-8) 




;<0KOfiiA4t, &*tr<8)X,li, <£#CxV§iftT&v>fc«;£5ll-Cl.>6. ztlbn 
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A 


B 


C 






0.600 


0.300 


0.200 




0.60J 


0.39S 


0.201 






O.SOO 


0.600 


0.000 




0.501 


0.599 


0.000 



B : Tot^uy-*/) -> 
c : <y!>JS-s> 



rv-4tv>ii. ifc, Bll.SliSll.lOfitffluTiE^&x-:* I- A<<8t8 
»l»30»-C*-jfc. 

11.4 ±»&fr&#rCDfSfll 



11.4 atMMMfWBB 22i 

bii.s as**? m. 




jr,= Se,+ u :/ = l.~.L (n-6) 

fc*>. ent h^T«-frBfi--»-i)it*^cfi. -9*9, Lmnmtx-i 



^ ar„ ar lf -, xi frrf,^C^$ii6ft*( JT t*«>6. 
X = ^» (11 .7) 

11.5 £fiS^tf£©:7D^5Zj<h3SgiiGfl 

► BASIC ?U 4? A. 

EU.7 iaWttVo/^A • tjji h 



10100 '*««****#**•••*»••*«**«»* 
ioie ' » « 
1020 *• Power Method « 
1030 '« « 
1040 '•***•*•«««•»»*•-««•»«•«*« 
1030 ' 

106O ' **• MAIN PROGRAM 
1070 ' 

1080 DEFJNT I-NlDEFDRL A-M.O-2 
1O90 CLS 3 

nee - 

111B lNPUT-No. of mxturn o -.ND 

1120 INPUT"No. a* Bamplino Points ° -.NS 
1130 INPUT H Data File nam' «= -.MMM* 



ii.5 &s.aa&£toTa73&tmnxm 

1140 IKPUT-Ho. of Eigenvalues - -,NE 

sisa ' 

1160 ID~«I>— 1 1 1&*1*>-1 : IE-KE-Ij J-0.K-O 
1178 ' 

1168 Din D(ID.I9>.CO1lD,X0> 
1 1 vo • 

1200 GOSUB «DIN 
1210 B0SUB «C0V 
1220 GOSUB *P0W 
1230 BOSUB •BRP 
1340 STOP 
1290 ' 

1260 ' *» DATA INPUT ** 
1278 - 

1280 *D1N t , , 

1290 Din Ml* (ID) 

1300 K-liFOT 1-0 TO ]D-liJ»INSTR<K.NMtt».-i-> 

1310 hff»II)BniD»(Krtl»,K,S*K)iK-OMiKeXT 

1320 MH» IID> «HI0* (KMM* ,K,LENtMHMt) — K+l ) 

1330 FOR 1-0 TO ID: OPEN "2i -♦WlOtI ) FOR INPUT AS 01 

1340 FOB O-O TO IGi INPUT* I ,D{ I . J) ■ NEXT 

1350 CLOSE*HNEXTtPRINT 

1360 RETURN 

1370 • 

1380 ' *• COVAftlANCE MATRIX 
1390 ' 
1400 *COV 

1410 FOR 1-0 TO I D t FOR J-l TO IDiCT-O 
1420 FOR K=0 TO ISiCT-CT+D ( I , K) »0( J.K) JNEXT 
1430 CO(l,J)-CT/NSiC0fJ,I>-CT/NStNEXT J.X 
1440 RETURN 
I43B ' 

1460 ' »* POWER ** 
147B " 
1480 *POW 

149B Din E(IE.ID) ,0<IE> ,VO<lD> .VNtID) ,EV(IE.IS) 
1SO0 EGP-.OOOBlitBP*200 
1310 HO»=0 
1S20 r-o 

1530 FOR 1=0 TO IDtVO(I)-ltNEXT 
1S40 IC-O 

ipso IC-IOl 

136© FOR 1-0 TO IDiW-OiFTJft J-O TO ID 
1570 W°M+CO(I.J)«V0tJ>tNEXTjVW<I)=W:NEXT 

1380 X-OiFDR I-O TO IDiH-VN(I) i IF ABS <M> >AB8 (X) THEN K-U 
1590 NEXT 

1600 FOR 1-0 TO IDrVNt I )-VN( 1 ) /XiNEXT 

1MB W*ABS(2/X-l)tIF W>£SP THEN 1660 

1620 FOR I**0 TO ID: IF ADS <VN 1 1 1 1 *ESP 1HEM l«4tt 

163Q fc*-ABB4VO(I)/VNO>-l> i IF U>EGP THEN 1660 

1640 NEXT 

163(3 GOTO 1690 

1660 IF IO=I6F THEN 1790 

1670 Z-XiFOft l-O TO ID* VO 41 > — VN < 1 1 i NEXT 
166B BCTO . 1550 

J690 X«O»F0R 1-0 FO lDiW-OtFOR J«0 TO ID 
1700 »-W+C0U.J1«VN(J)iKEXTiX«X+VMU)«M:NEXT 
1 71 B WOt FOR 1-0 TO 1DiW-W*VNI1)«VN<I»iNEXT 
172G X-X/UiG(NO)-X 



173© FOB 1-0 TO lDiE<NO.I)-VN(ltiKEXT 
1740 IF W0>«=IE THEN 1820 

J 730 W-ClFOR 1-0 TO I Di < I ) -VM < I ) i-KEXT 

I7U FOR 1-0 TO IDt Z— VNU ) *X/« 

1770 FOR TO IDsCOd . JI-CO(I, J)«VN( J) *ZiNCXT J. I 

1760 NO»ttO*liGOTO 1520 

1790 PRINT" ** CALCULATION DOSE NOT CONVER6ENCE 
1BOQ 6 (NO) -C: FOR I-Q TO ID«E(ND ( Il-OiNEXT 
1810 ' 

1830 FOR 1-0 TO IEiGMl-e.FUR Joe TO ISiSN2— O 
1630 FOR K=C TO lDiSM2=SM2*-Etl ,K>*D(K. JI iKEXT 
1640 EV<I .J)-6M2iBHl-Sm*6M2»SH2«NEXT 
18S0 SH1«SCR<SM1) 

FOR J-0 TO 15iEV<I,J>-£Vli,J>/BmiNEXT J.I 
RETURN 

•« GRAPHIC *« 



1870 



1900 
1910 
1920 
1930' 
1940 
1950 
I960 
1770 
1980 
1990 
2080 
2010 
2030 
2030 
2040 
2050 
2060 
2O70 
2080 
2090 
210* 
2110 
2120 
2130 
2140 
21B0 
2160 
2170 
2160 
2190 
2300 
2210 
2220 
2230 
224B 
2230 



CONSOLE ..O.OtSCREEN SfCLS 1 
PRIMT>PRlNTtPRlMT a Ci9«nv*luM a 
FOR 1=0 TO IE* PRINT 6UHNEXT 
DIN X0U>,Y0tl) 
XS-253I XD»XS/ 1 Bi YS- 1 50 

XO (0) -10» XO< 1 1-3101 YQ <0> «200i YO< 1 > -375 

LINE<X0(1) .VO(0> )-(XO(l>+Jl9,YOl0)-YS> ,7 ,B 
LINE<X0<1> ,YOUH-UOU)+XB,YO(l)-YB> ,7,B 
LINE CX0< I > . YO ( l>-YS/2) - CXO ( 1 > *X8 . YO < 1 > -YS/2> ,7 

AHAX«OlFOR 1=0 TO ID: FOR K-0 TO IS 
IF 0<I,K)>AMAX THEN ANAX«0(I.K> 
NEXT K.IlSC«YS*.9/AMAX 

FOR 1-0 TO IDlX-XOU>iY-VO<0»-6C»D<I.O) 

RBET(X,Y) ,4fFQB J=C TO IS 

X«il«XCH-X0 (Di Y-YO(O) -SOD « 1 . J » 

LINE — IX . Y) ,4i NEXT J.l 

BEET UFOR 1*1 TO 236 1 NEXT I BEEP <3 

PRINTS PRINT* INPUT "No. of Eigenvectors - -|N£ 
lE-NE-t 

AHAX=0rFOR !«8 TO IE t FOR K«0 TO IS 

IF ABB <EV ( I «K> ) >AHAX THEN AHAX-ABS(EV U ,K> ) 

NEXT K.lt8C»YS*.9/AHAX/2 

FOR 1-0 TO IEiX-XO<l>iY-YO<l)-Y6/2-eC«EV<I.0> 

P6£T<X,Y) ,4. FOR J°0 TO 18 

X-J»XD*XOtl) iY-Y011)-Y3/2-SC»EV(I ,Jt 

LINE -<X,Y),4aNEXT J, I 

BEEP llFOR I"! TO 2501 NEXT* BEEP 0 

RETURN 



► 9t*X0l 
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1200 948 1200 

» B (cnr 1 ) 

*£fefl4. COTitiO. 151SW Jf •> u >KHIiZ-3W»tt»(o-, m-, />-w8t£<*:)* , S 
al,T^>fcCt*•l3l>•■3^S■ , . 
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(3) M27"-;?e*slt4«erli. JbS6MS#t*9, e$fcli»dtC*4C: fc. 
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